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In supersymmetric (SUSY) models [H, the spin- 1/2 
quarks and leptons have spin-0 quark and lepton part- 
ners. Experimental data suggest that the superpartners 
of the first and second generation are heavier than those 
of the standard model (SM) particles, while the mass of 
the lighter scalar top quark (stop or t{) is weakly con- 
strained and can be below that of the top quark . This 
is due to the mixing between the left and right handed 
interaction eigenstates which is a function of the large 
Yukawa coupling of the top quark. At the Fermilab Teva- 
tron stop quarks and antiquarks can be produced in pairs 
in strong interactions (gg/qq — > t\ti). A single stop could 
also be produced at the Tevatron, e.g., via bg — > t\T [3j; 
however, unlike pair production, this process requires an 
i?-parity (R p ) violating vertex. In regions of parameter 
space not excluded by data, R p violating (ftp) couplings 
are small [H, making single stop production negligible 
compared to pair production. Stop quarks can decay 
into lighter SUSY and SM particles if R p is conserved 
or into ordinary quarks and/or leptons if R p is violated. 
Within the framework of ftp SUSY [4] , theoretical studies 
indicate that the dominant decay mode for the light stop 
is the lepton number violating decay t\ — > rb for a wide 
range of SUSY model parameters, including the region 
allowed by neutrino oscillation data 

Leptoquarks appear in various SM extensions Q. 
Charge \Q\ = 2/3 and \Q\ =4/3 third generation scalar 
leptoquarks (LQ3) are expected to decay into r and 
b with B(LQs — > rb) = 1 for all LQ3 states when 
m(LQ 3 ) < m(t). The next-to-leading order (NLO) cross 
section for LQ^LQ 3 production is very close to the titi 
production cross section cr(tiii) since diagrams with vir- 
tual gluino exchange are strongly suppressed with the 
existing limits on gluino mass [7j. Thus, the limits ob- 
tained for ftp stop should be fully applicable to the LQ 3 
case. 

In this Letter we describe a search for titi — > T + r~bb 
with the upgraded CDF detector (CDF II) || and set 
an upper limit on cr(titi) x /3 2 , neglecting additional de- 
cay modes that may pass selections of this analysis when 
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(3 = B(ti — > rb) < 1. We look for a final state with either 
an electron or muon from the decay r — > lvgy r (£ = e 
or /i), a hadronically decaying tau Th, missing transverse 
energy f}p [t| from the neutrinos, and two or more jets. 
We have studied the addition of a requirement that the 
jets are consistent with originating from the hadroniza- 
tion of a b quark but found that the increase in purity is 
outweighed by the loss in signal acceptance. Therefore, 
we make no such specific requirement. This analysis uses 
approximately three times more data at a higher y/s than 
the previous CDF result that set a 95% C.L. limit 
of m(ti) > 122 GeV/c 2 . The increased yfs is expected 
to give a substantial increase in the t\ti production rate, 
e.g., - 35% for m(t x ) = 155 GeV/c 2 . 

CDF II features several main subsystems critical to 
this analysis. The charged particle tracking system con- 
sists of multi-layer silicon detectors and an open-cell 
cylindrical drift chamber enclosed in a 1.4 T supercon- 
ducting magnet. At |ry| < 1 [9j charged particle trajec- 
tories traverse all chamber layers, while at larger |ry| the 
chamber coverage is reduced progressively. The calorime- 
ter system is organized into electromagnetic and hadronic 
sections segmented in a projective tower geometry and 
covers \r}\ < 3.6. A set of strip and wire chambers is 
located within the central electromagnetic calorimeter at 
approximately the depth of shower maximum and aids in 
reconstructing electrons, and photons for |ry| < 1.1. The 
muon detection system is located outside of the calorime- 
ter and covers \rj\ < 1.0. 

The analysis begins with a data sample collected by 
inclusive lepton-plus-track triggers [ll|. These triggers 
select events containing an electron (muon) candidate 
with Et > 8 GeV (pr > 8 GeV/c) and a second track, 
which is required to be consistent with originating from 
a tau decay by demanding that there be no other nearby 
tracks with px > 1.5 GeV/c between the cones of 0.175 
and 0.524 radians around the track. The integrated lu- 
minosity of the data sample is 322 ± 19 pb _1 [l2|. 

From the trigger sample we select events offline by 
identifying at least one lepton with pt £ > 10 GeV/c and 
at least one Th candidate in < 1. The details of the 
Th identification algorithm can be found in Refs. [l3l.[l4j]. 
We require p T T > 15 GeV/c. Jets are reconstructed using 
a fixed-cone algorithm with AR = \J A77 2 + A0 2 = 0.4 
within < 2.4. 

The dominant SM backgrounds are vector boson pro- 
duction, QCD, and ti production. In QCD multijct 
events, for example, scmileptonic b quark decays or 7 con- 
versions can be misidentified as lepton candidates, and 
narrow jets can be misidentified as Th candidates. Wc 
require the sum of the px of the tracks within AR < 0.4 
around the lepton candidate (Itrk) be less than 2 GeV/c, 
and no jet with E T > 15 GeV within 0.3 < AR < 0.8 
around the lepton. Further, we reject events where the 
muon or electron candidate is consistent with a cosmic 
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ray muon or 7 conversion electron (see Ref. [13] for de- 
tails). We veto events where the invariant mass of the 
primary electron (muon) and a reconstructed electron 
(muon) candidate,which is required to pass loose iden- 
tification criteria [l3|, is 76 < mu < 106 GeV/c 2 . We 
also reject events with 76 < m eT < 106 GeV/c 2 and az- 
imuthal separation of |A0 er | > 2.9 rad. For the muon 
channel we do not apply a similar requirement, as the 
probability for a muon to be reconstructed as a Th is neg- 
ligible. To suppress further QCD and Z — > tt events [lOj, 
we require S T = Wt\ + \Pr h \ + \Pr\l c > 110 GeV/c. 
We define six regions in the mT(£,fr) = 

^j2"p l T fr(\ — cos A^^) versus iVj e t plane, and denote 

them as Aj (B 3 -) for m T < 35 GeV/c 2 (m T > 35 GeV/c 2 ) 
and j = 0, 1 or 2 for iV jct = 0, 1 or > 2. We count 
into iVjet the jet candidates that have Et > 20 GeV 
and are separated from any of e, /i or Th by AR > 0.8. 
The minimal values of St and jet Et are optimized for 
maximum significance in the A2 region for 140 — 160 
GeV/c 2 ti's. The m T < 35 GeV/c 2 cut effectively 
separates signal from W + jet and it backgrounds. The 
^Vjet > 2 requirement strongly suppresses the Drell-Yan 
contribution. The data in region A2 were not examined 
until the analysis procedure was finalized. The regions 
with iVj e t = or 1 contain mostly background and were 
used mainly as control samples for validation. Region 
B2 has an appreciable signal acceptance (~ 40% of 
that in region A2) but substantially higher background 
expectation. For statistical interpretation of the data, 
we developed a likelihood method that, in addition to 
our primary signal region A2, utilizes side-band regions 
Ao, Bo, and B2, which are used to perform data-driven 
W + jet background estimations and to improve the 
sensitivity of the analysis. 

The total event acceptance is a = 6mc 1 e trig- Here 
£mc is the product of geometrical and kinematical accep- 
tances, efficiencies to identify lepton and Th candidates, 
including isolation requirements, and the efficiency for 
the all remainin g cu ts. We use pythia Monte Carlo 
(MC) generator and the GEANT3-based [ll] CDF 
II detector simulation to calculate smc- Our nomi- 
nal choice for the parton distribution functions (PDFs) 
is CTEQ5L [13] with the renormalization scale Q = 
^/m(ti) 2 +pT(t\) 2 . The trigger efficiency e tr i g is mea- 
sured using data [l4j |. In region A2 a increases nearly 
linearly from about 0.6% at m(ti) = 100 GeV/c 2 to 2.7% 
at 170 GeV/c 2 and is similar for both electron and muon 
channels. 

The combined systematic uncertainty on a decreases 
almost linearly from 11% for m(h)= 100 GeV/c 2 to 7.2% 
for 170 GeV/c 2 and is similar in both channels. The 
largest contribution comes from the PDF systematic un- 
certainty, which is estimated using the uncertainty sets 
of CTEQ6.1M PDFs and the technique described in 
Ref. . For a 150 GeV/c 2 stop this uncertainty on a is 



4.0%. The uncertainty due to an imperfect knowledge of 
the jet energy scale, determined by varying the scale by 
±lcr, is 2.9%. The uncertainty due to the amount of ini- 
tial and final state radiation is found to be 2.5%. Other 
sources of systematic uncertainty include the uncertain- 
ties in lepton and Th identification and isolation, and fr 
resolution, and amount to a 5.1% relative contribution. 
The uncertainty on the integrated luminosity is 6%. 

The SM backgrounds come from two sources: (i) events 
with a true ir^ pair from Z/^*(~+ Tr)+jets, it and dibo- 
son (WW, WZ, ZZ) production; and (ii) events where 
lepton or Th candidates do not originate from a true lep- 
ton or tau but from the jets in W + jet, Z/j*(-^ l£)+]ets 
and QCD multijet events. We first estimate the back- 
ground from SM processes excluding the W + jet con- 
tribution. Drell-Yan, it, and WW production are es- 
timated using pythia and the CDF II detector simu- 
lation. For Drell-Yan backgrounds we use scale factors 
that improve the agreement between the prediction for 
the yield of these events in MC simulation and the yield 
observed in data. The QCD multijet contribution is es- 
timated by extrapolating the number of observed events 
in data for events with non-isolated leptons, defined by 
2 < Itrk < 10 GeV/c, into the class of events with an iso- 
lated lepton, defined by I trk < 2 GeV/c 0. The NLO 
cross sections of 6.7 ± 0.7 pb [lj| and 13.5 ± 0.5 pb [2(| 
for it and WW production, respectively, are used. The 
contributions from WZ and ZZ are found to be negligi- 
ble. 

The pythia MC simulation does not accurately pre- 
dict the absolute rate of the W + jet background con- 
tribution (N w ) in this analysis. Therefore, to esti- 
mate N w in each region, we use the differences be- 
tween the data and all other backgrounds plus signal 
in regions A2, B2, Ao, and Bo and the assumption that 
U = [N W {A 2 )/N W {B 2 )] ■ [N W {B )/N W (A )] ~ 1. The 
ratios in 1Z are determined by kinematics of the W + jet 
events at fixed Nj t and are well modeled in MC. Based 
on MC predictions and cross checks with data vs MC 
comparisons, we conclude that 1Z = 1.0 ± 0.5 is a conser- 
vative assumption. We define a likelihood function using 
Poisson statistics as a function of a (hi 1) and N w . The 
input parameters to the likelihood are the numbers of 
observed and expected events in each of the four regions. 
The number of expected events in region i is given by 
Ni = a(hh) ■ B(tt -> ^T h ) • JCdt- ai + N? G + NW, 
where the branching ratio B(tt — > It\^) ~ 0.23, Nf G in- 
cludes all SM backgrounds except W + jet events, and 
a,; is the total event acceptance for signal in region i. 
The ratio 1Z = 1.0 ± 0.5 and N w in regions A , B , and 
B2 are nuisance parameters with flat prior distributions. 
The large uncertainty on TZ does not affect our final re- 
sults because N W (A 2 ) is expected to be small. We use 
this two-dimensional likelihood to estimate N w for each 
region and to calculate upper limits on o(hh) x 1 . 

In Table U we show the number of events observed 
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TABLE I: Number of events observed in data, iV b 8 , along 
with the expected number of SM background events. The 
W + jet contributions are shown separately as they are esti- 
mated using the observed number of events in the data, the 
SM prediction excluding the W + jet contribution, and a pos- 
sible stop quark contribution. 
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FIG. 2: 95% C.L. limit curves for the titi production cross 
section times the branching ratio t\ — > br squared for the 
cases when the theoretical uncertainty on the cross section is 
(dashed line) and is not (dotted line) considered in the limit 
calculation (see text for details). The previous constraint 
m{LQ 3 ) > 99 GeV/c 2 [HI is also shown. 



TABLE II: 95% C.L. upper limit on a{t x ix) x (3 2 (in pb) 
as a function of m(ti) for the cases when uncertainty on the 
theoretical cross section is considered {a^^ h uncert x /3 2 ) and 
is not considered {p n ^o°u.ncert x /3 2 ), where j3 = B(ti — > rb). 



m(ti) (GeV/c 2 ) 100 110 120 130 140 150 160 170" 

°w%k uncert * P 2 (pb) 4.73 3.37 2.50 1.99 1.61 1.38 1.26 1.14 
ofo uncert x P (p b ) 4.48 3.11 2.27 1.81 1.47 1.26 1.16 1.04 



FIG. 1: Distribution of the number of jets (E T > 20 GeV) 
for events with irt^It) < 35 GeV/c 2 (regions Ao, Ai, and 
A2) compared to the expectations from SM processes and 
prediction for titi (m(ti) = 150 GeV/c 2 ) signal. 



in the data along with the SM expectation. In Fig. Q] 
we present the iVj e t distribution for events with ttit < 
35 GeV/c 2 (regions A , Ai, and A 2 ). 

Two events are found in region A2, consistent with the 
SM prediction. We use the likelihood function to ob- 
tain a probability of such an observation given a specific 
signal cross section, and calculate a 95% C.L. limit on 
<r{t\t\) x 1 . The electron and muon channels are treated 
as two separate measurements, taking into account the 
correlations among the systematic uncertainties. 

Figure [2] shows the 95% C.L. limit curves for a(tit\) x 
/3 2 as a function of m(ti), with the numerical results given 
111 Table HH The dotted curve is our experimental result, 
compared to the NLO cross section (solid line) obtained 
using PROSPINO version 2 [23] with our nominal choice 



of CTEQ6.1M PDFs [l3,_and Q. The theoretical un- 
certainty of ±18% on a{t\ii) is due to the choice of Q 
(varying the scale from its nominal value by a factor of 
two or a half) and PDFs. Taking this uncertainty into 
consideration, the limits are re-evaluated and are shown 
in Fig. [2] using a dashed line. The corresponding mass 
limits for the first and second cases are 156 GeV/c 2 (com- 
pared to 122 GeV/c 2 0) and 153 GeV/c 2 respectively. 

In conclusion, we have searched for t\t\ production 
in the final state of a lepton (e or fi), a hadronically 
decaying tau, and two jets using 322 pb" 1 of pp collision 
data at y/s = 1.96 TeV. We observed no excess of events 
in data over the number of expected events from SM 
processes. In an ^tp SUSY scenario, we set a 95% C.L. 
lower limit on the t\ mass of 153 GeV/c 2 taking into 
account the theoretical uncertainties on the NLO cross 
section and assuming B(t\ — > rb) = 1. These results are 
also applicable to LQ3 pair production. 
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